Critical dynamics of singlet excitations in a frustrated spin system 
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We construct and analyze a frustrated quantum spin model with plaquette order, in which the 
low-energy dynamics is controlled by spin singlets. At a critical value of frustration the singlet 
spectrum becomes gapless, indicating a quantum transition to a phase with dimer order. The 
magnetic susceptibility has an activated form throughout the phase diagram, whereas the specific 
heat exhibits a rich structure and a power law dependence on temperature at the quantum critical 
point. This generic behavior can be relevant to quantum antiferromagnets on Kagome and pyrochlore 
lattices where (almost) all low-energy excitations are spin singlets, as well as to the CaV409 lattice 
and the strongly frustrated antiferromagnets Cu2Te205(Cl,Br)2 and Li2VO(Si,Ge)04. 
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It was argued by Haldane and Chakravarty, Halperin 
and Nelson that a quantum phase transition between a 
Neel state and a magnetically disordered phase in a two- 
dimensional (2D) quantum antiferromagnet is described 
by the 2-1-1 dimensional non-linear 0(3) cr-model. More 
recently the problem was studied in detail using the 1 /N 
expansion {N = 3 is the number of components of the 
order parameter), numerical methods |^ as well as the 
Brueckner method In this description the lowest 

quasiparticle excitation is a spin triplet. There are cases 
when a low-energy singlet excitation also appears, how- 
ever it turns out to be irrelevant to the critical dynamics 
because of its vanishing spectral weight at the quantum 
critical point Q. 

The above picture based on S=l excitations describes 
a wide class of quantum phase transitions j7|. On the 
other hand in Cr-based S=3/2 Kagome materials it was 
found that at low temperatures (T < 3 — 5 K) the specific 
heat C oa and is practically independent of magnetic 
field up to 12 T If we try to understand these data 
in terms of quasiparticles, the temperature behavior in- 
dicates a linear dependence of the excitation energy on 
momentum, and the insensitivity to field suggests that 
the quasiparticle spin is zero. Finite cluster numerical 
simulations for the S=l/2 Kagome model indeed show 
that the lowest excitations are spin singlets with very 
small or zero gap Similar behavior with many sin- 
glet states inside the triplet gap was also found in the 
pyrochlore antiferromagnet and the CaV409 lattice 
PI . In the Kagome case the singlet quasiparticle picture 
is far from being certain. It is unclear how to explain in 
this scenario the rather large zero temperature magnetic 
susceptibility observed in the Kagome materials |^ . It is 
also possible that the material is a spin glass and hence a 
description in terms of quasiparticles would not be ade- 
quate. Nevertheless, it is very important and interesting 
to analyze the quasiparticle scenario. From this point of 
view the Kagome magnets are critical (or close to crit- 
ical) systems, and they certainly can not be described 
by the 0(3) cr-model. In the present work we consider a 



two-dimensional spin system which exhibits a non-trivial 
singlet dynamics leading to a quantum critical point con- 
trolled by singlet excitations. 

Critical singlet dynamics can naturally appear in a 2D 
quantum spin model with a zero temperature quantum 
phase transition separating two singlet ground states, 
both of which are magnetically disordered but have dif- 
ferent discrete lattice symmetries. We show that a spin 
model constructed as an array of weakly coupled frus- 
trated plaquettes has only singlet low-energy degrees of 
freedom in a certain range of parameters. Many the- 
oretical works have been devoted to the conventional 
order-disorder quantum phase transition in similar mod- 
els, mainly in relation to the CaV409 compound [|2|-0. 
Different types of singlet dimer states have also been dis- 
cussed for the quantum dimer model! 15[ , the Heisenberg 
model on a 3D pyrochlore lattice |ll], and the general- 
ized 2D Ji — J2 model However, low-energy critical 
singlet dynamics for these systems has not been stud- 
ied until now. Very recent numerical work has shown, 
in the context of the 2D Ji — J2 model, that a quantum 
transition governed by singlet dynamics takes place for 
strong frustration J2 ~ Ji/2 jlj]. It was also recently 
argued that this regime is experimentally accessible in 
Cu2Te205(Cl,Br)2 and Li2VO(Si,Ge)04 0, and mea- 
surements of the excitation spectrum of these materials, 
such as Raman spectroscopy, are now actively being pur- 
sued [ pO[ . Let us mention that dimer states appear also 
in various ID spin models (e.g. in the frustrated J1 — J2 
Heisenberg chain). Their critical behavior is understood 
in terms of deconfined spin-1/2 excitations or spinous, 
which do not generally exist in higher dimensions. A ID 
system with features similar to the generic critical be- 
havior we analyze in this work is an effective ID Kagome 
lattice model [|2l|] with gapless singlet and gapped triplet 
excitations. However, unlike the behavior we focus on, 
this model does not exhibit a quantum phase transition 
and breaking of lattice symmetry. Thus our future dis- 
cussion is relevant to models in D > 1, and a 2D case is 
analyzed as a generic example. 
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Consider a cluster of four S=l/2 spins (plaquette, see 
Fig. 1(a)): 

Ho - Ji^S,-S, + J2^S.-S„ (1) 

where {ij) and (ij) denote, respectively, the side and di- 
agonal bonds. We choose the diagonal couphng J2 to 
be close to the side one: J2 « Ji = J > 0, and de- 
fine a — Ji ~ J2 ^ J- In this regime the spectrum 
of a plaquette is the following: two close singlet states 
which we denote by \sa) and \sb), with energy differ- 
ence 63 — £A = 2a, three almost degenerate triplet states 
with excitation energy J above the singlets, and one 
S=2 state with energy 3J. The two singlet wave func- 
tions are expressed as: \sa) = ^{[l,2][3,4]-h[2,3][4, 1]} 
and \sb) = {[1,2] [3, 4] - [2, 3] [4,1]}, where [i,j] denotes 
a singlet formed by the nearest-neighbor spins i and j 
(Fig. 1(a)). Both singlets \sa) and \sb) are invariant un- 
der a four- fold rotation of a plaquette, whereas the colum- 
nar dimcr states, e.g. [1,2] [3, 4], break this symmetry. 
Consider now two decoupled plaquettes which for a — 
have a 2 X 2 degenerate ground state. The degeneracy 
is lifted if a weak interaction between the plaquettes is 
switched on selecting one of the columnar dimer states. 
Such a degeneracy lifting is different from the effect of a 
non-zero a and competition between the two leads to a 
nontrivial singlet dynamics in an array of weakly coupled 
plaquettes. 

To be specific consider a square array with antifer- 
romagnetic couplings ji between nearest-neighbor spins 
and j2 between next-nearest-neighbor (diagonal) spins 
from different plaquettes, see Fig. 1(b). For ji — Ji and 
j2 = J2 it would be equivalent to the translationally- 
invariant 2D Ji — J2 model. However now we are inter- 
ested in the weak-coupling limit: ji,j2 J2 ~ Ji = J ■ 
The low-energy singlet sector of the Hilbert space has a 
natural pseudospin representation in terms of the states 
IT) = \sa) and 1 1) = \sb)- The total Heisenberg spin 
Hamiltonian of the array of coupled plaquettes is mapped 
on the following pseudospin Hamiltonian to lowest order 
in the small parameters a and ji.2/ J'- 

nj (-) \ ^ QZ 1 1 QX QX I 6^*^^'^^^ / nZ QX i QX QZ \ 

^nhJ2s: . (2) 

i 

The energy scale f2 and the effective dimensionless "mag- 
netic field" h are obtained in second-order perturbation 
theory: 

" = ^(^1 + ^2), = ^[2« + ifi + jI - 6jij2)/6J]. 

(3) 

We have also defined Q = (0, tt) and dropped a constant 
term. The first sum in Eq.(p|) is over nearest-neighbor 



plaquettes since bonds connecting second neighbors con- 
tribute to a constant term only and thus do not change 
the singlet dynamics. We stress that (^) is an exact map- 
ping of the original Heisenberg model in the low-energy 
sector (excitation energy w ~ ri,a ^ J). The Hamilto- 
nian Ti. describes an anisotropic ferromagnet in an exter- 
nal field. The interaction between pseudospins on adja- 
cent sites is purely Ising, which can be seen by rotating 
the spin axes. However the Ising axis is staggered in the 
x-z plane deviating by angle ±7r/6 from the x-direction 
for horizontal (vertical) pairs. 

The form of the Hamiltonian Eq.(2) is unaffected by 
asymmetry in the coupling constants ji and j2, which 
would only change the position of the reference zero- 
field level. In fact, the pseudospin Hamiltonian remains 
the same with parameters given by: Q = jf/2J and 
h = [2a + {ji — 3jij2)/6J]/ri, even if we switch off one 
of the diagonal j2-bonds between nearest-neighbor pla- 
quettes. In this case our model resembles the frustrated 
^-depleted square lattice of CaV409 [|l^-[l^, which has 
next-nearest-neighbor couplings comparable or even ex- 
ceeding the nearest-neighbor exchange. Thus our results 
apply to the singlet ground states of this magnet as well. 

In zero "magnetic field" = the ground state of 
7i has broken Ising symmetry with all spins parallel or 
antiparallel to the x direction. In the language of spin 
singlets breaking of Z2 symmetry corresponds to sponta- 
neous dimerization in one of the two columnar patterns. 
The rotational symmetry of the square lattice is obvi- 
ously broken. The "magnetic field" tends to orient all 
spins along the z axis, and depending on the sign of h 
favors either \sa) or \sb)- Hence there are two critical 
fields, a positive and a negative one, for transitions into 
states with a restored Z2 symmetry. Since Eq.(2) is in- 
variant under h —h, we consider the case h > only. 

To study in more detail the properties of the symmet- 
ric phase near the critical point we map the pseudospin 
Hamiltonian (H) onto a hard-core boson model. Positive 
and large h corresponds to Ji > J2. In this case we 
choose the bare ground state as |0) = JJ- \sA)i- A boson 
creation operator on site i is defined by \sB)i — bl\sA)i- 
In terms of pseudospins we have Sf ~ \ — h\hi, — b]. 
The resulting boson Hamiltonian is: 

Hb =eY,blb. + tY.iblb, + bib] + h.c.) + (4) 

i (ij) 

g Hb\b]b. + b\b% -f h.c.) + y ^ b\b]b,b,. 

In the (jf-term the sign -|- corresponds to a horizontal link 
and the sign — to a vertical one. The parameters are: 

e = n{\ + h), t=-f , 5 = ^ = -f . (5) 

Considered in combination with the hard-core constraint 
{b\Y = 0, the Hamiltonian (|^) is equivalent to the pseu- 
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dospin Hamiltonian and hence it is an exact mapping 
of the original spin problem in its low-energy sector. 

The bosonic form of the effective Hamiltonian is con- 
venient for the analysis of the low density (disordered) 
phase. Diagonalization of Eq.(4) in the quadratic approx- 
imation gives the following excitation spectrum: cjk = 
+ 4t7k)^ - (4t7k)^, where 7k = ^(cosfc^; -I- cos ky). 
We set the inter-plaquette lattice spacing to one. The 
excitation gap A = a;k=o vanishes at the critical point 
/ic = |. The other critical point is at —he- For h > he 
the model is in the A-type plaquette phase, for h < —he 
it is in the B-type plaquette phase, and in between 
—he<h<he there is a phase with a boson condensate 
at k = 0. The condensate is doubly degenerate because 
the boson field can have both signs. 

The zero-point quantum fluctuations, which exist in 
the disordered (plaquette) phases, change slightly the 
critical field he- The simplest way to account for the 
correlation effects is to use the Brueckner technique de- 
veloped in Ref. Q . The cubic and quartic interactions in 
(^ are treated in the one-loop approximation, and the 
hard-core constraint is enforced via an infinite repulsion 
term Uj^i^l^i^i^ij U oo, added to the Hamiltonian 
(^ . This interaction is taken into account in the Brueck- 
ner (low-density gas) approximation. The singlet gap 
obtained from the self-consistent solution of the result- 
ing Dyson equations is plotted in Fig. 2. At the critical 
point the small parameter which justifies the Brueckner 
approximation is the singlet density (blbi) « 0.04. The 
renormalized critical value he ~ 0.65 is only slightly be- 
low the result for non- interacting magnons. The reason 
is the strong compensation of the hard-core corrections 
by the one-loop diagrams arising from the cubic terms in 
Eq.(4). The dependence of the gap on h near the tran- 
sition point agrees with the critical index v = 0.63, ex- 
pected for the 0{N = 1) (t model and is quite differ- 
ent from the value ly — 0.5 for non-interacting magnons. 

Next we analyze the singlet excitation spectrum in the 
ordered dimer phase, h < he- We return to the pseu- 
dospin representation ^ and use the analogy between 
the evolution of the dimer phase and the spin reorienta- 
tion process in an external magnetic field. The Holstein- 
Primakoff transformation is applied to the pseudospins in 
the rotating coordinate frame, which to lowest order coin- 
cides with the hard-core boson representation. At h = 
all spins point along the x-axis. For finite h the spins 
tilt towards the field direction at an angle sm6 — h/he- 
Keeping only quadratic terms we find the following "clas- 
sical" spectrum of singlet excitations: 

= r!^l - 7k(l - I cos2 e) + ^ sin207k+Q . (6) 

In accordance with the broken discrete symmetry the 
singlet excitations have a finite gap A = Wk=o — 
i7-y/2(l — /h1)/i, plotted in Fig. 2. The gap vanishes at 
the critical point. Proceeding further with the spin anal- 



ogy we have also calculated the "spin reduction" for the 
ferromagnet (H). This parameter shows how far is the 
real ground state wave-function from the approximate 
mean-field ansatz. At /i = we find {S) ~ 0.498, mean- 
ing that the effect of quantum fluctuations is extremely 
small and the linear spin-wave theory is well justified. 
At the transition point h = he the zero-point oscillations 
are larger (S) « 0.44 but still small enough to justify our 
approximations leading to Eq.(||). 

The broken Z2 symmetry in the ground state at \h\ < 
he leads to a finite temperature transition. The whole 
phase diagram of the system of coupled plaquettes in the 
h-T plane is shown schematically in Fig. 3. The crit- 
ical temperature can be estimated as Te{h = 0) ~ 
(1.14r2 for a pure Ising case). The phase transition be- 
tween the ordered dimer and the plaquette states be- 
longs to the 2D Ising universality class, and we there- 
fore expect a logarithmic singularity in the specific heat: 
C ~ In \T—Te{h)\. Below Te the excitations acquire a gap 
and C goes exponentially at low temperatures. The same 
activated dependence holds also for the symmetric phases 
at \h\ > he- However, when we approach the quantum 
critical points at h = ±he the gap becomes smaller and 
the low-T behavior of the specific heat changes to the 
quantum critical law Cy = jT'^. The universal prefac- 
tor 7 can be calculated using the Brueckner technique 
g. The resuh is 7 = « per plaquette, 

where c = il/2 is the singlet excitation velocity at the 
critical point (we set fcs = h — 1). This value coincides 
with the large-N mean-field result Q. Notice that, un- 
like Cv, the magnetic susceptibility is always activated 
X ^ exp(— At/r),r 0, since it is governed by the 
triplet gap At ~ J. 

The variation of Cy as a function of temperature is 
schematically presented in Fig. 4 for the different parts 
of the phase diagram. In addition to the singlet con- 
tribution, whose form was discussed above, we have also 
shown the peaks expected to arise from the higher energy 
triplet and quintiplet states. We remind the reader that 

and J represent two distinctly different energy scales 
since according to our weak inter-plaquette coupling as- 
sumption J. We have estimated that the singlets 
contribute about 25% to the entropy S = J {Cv/T)dT 
(and 61% and 14% are taken by the S=l and S=2 states, 
respectively). The possibility of rich behavior at low T 
has been debated for some time for the Kagome antifer- 
romagnet and a sharp structure is believed to exist due to 
low-energy singlets . Notice that in our model a very 
peculiar behavior with a logarithmic singularity sets in at 
h < he which crosses over to the quantum critical regime 
at h = he- Such a sharp feature is generic for systems 
near a transition between two singlet ground states. We 
also note that near the quantum critical point h w he, Te 
is expected to be small (and ultimately vanish at h — he) . 
Consequently the logarithmic singularity at T = Tc ^ 
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and the singlet peak at T ~ should be clearly sepa- 
rated (this has been assumed in Fig. 4). Alternatively, for 
h 1 one has Tc ~ and therefore the two contributions 
should merge. A numerically reliable calculation of Cy 
in the temperature region T < is a separate problem, 
however for higher temperatures (of order J) a quanti- 
tative description can be easily achieved, and the result 
is plotted in Fig. 5 (for the specific value of O = O.IJ). 
The contributions of the triplet and quintiplet states have 
merged into a broad peak around T « J/3, while the 
singlet peak is sharper and centered at T « $7/2 (in this 
region our calculation is expected to be qualitatively cor- 
rect). 

In conclusion, we have analyzed a quantum spin model 
with purely singlet low-energy dynamics. There is a 
quantum phase transition in the model separating a dis- 
ordered plaquette phase and a columnar dimer phase. 
Even though the inter-plaquette interaction was assumed 
to be weak, we expect our results to hold also for stronger 
interactions as long as there are no other instabilities, and 
to be applicable, e.g. to CaV40g. The broken symme- 
tries in the singlet ground state of CaV409 were con- 
sidered previously in Ref. JTst in the framework of the 
quantum dimer model, and the possibility of Ising tran- 
sitions between spin-Peierls and other disordered phases 
was discussed [ p4[ . Furthermore, recent numerical stud- 
ies of the square-lattice Ji — J2 model have indeed 
found a quantum transition near J2 ~ Ji/2 of the type 
discussed in the present work, even though a formally 
small expansion parameter (justifying the separation of 
the singlet and triplet dynamics) is not available. 

We have found that while the magnetic susceptibility 
is always activated, the specific heat behavior is very rich 
and changes substantially in the different regimes. The 
model was inspired by experimental and numerical data 
for Kagome systems and also shares common structure 
with pyrochlore antiferromagnets, the CaV40g lattice, 
and strongly frustrated square lattice antiferromagnets. 
The novel quantum critical behavior associated with sin- 
glet criticality (or proximity to such a critical point) dis- 
cussed in this work can be relevant to a wide class of 
disordered quantum spin systems, and could be detected 
by measurements of the specific heat as well as the low-T 
excitation spectrum in the S=0 sector via Raman spec- 
troscopy. 
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FIG. 1. (a) A single frustrated plaquette. (b) An array of 
interacting plaquettes. 
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FIG. 2. Singlet gap in the vicinity of the quantum critical 
point separating the plaquette and dimerized phases. 

FIG. 3. Schematic phase diagram of the model. The solid 
line represents the 2D Ising critical line, and the dots are the 
two quantum critical points belonging to the 0(1) (equiva- 
lently 3D Ising) universality class. The vertical dashed line 
represents the crossover boundary between the two plaquette 
phases. 

FIG. 4. Schematic temperature dependence of the specific 
heat throughout the plaquette phase, h > he, down to the 
quantum critical point h = ±hc (top), as well as for h < he, 
across the 2D Ising critical line (bottom). 

FIG. 5. Calculated specific heat for O = J/10, showing the 
low-temperature singlet contribution at T « n/2 as well as 
the broader peak on the scale of J arising from the S = 1,2 
states. 
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